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ABSTRACT 

We present the redshift distribution of the SHADES galaxy population based on the 
rest-frame radio-mm-FIR colours of 120 robustly detected 850/im sources in the Lock- 
man Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift 
distribution derived from the full SED information is shown to be narrower than that 
determined from the radio-submm spectral index, as more photometric bands con- 
tribute to a higher redshift accuracy. The redshift distribution of sources derived from 
at least two photometric bands peaks at z w 2.4 and has a near-Gaussian distribu- 
tion, with 50 per cent (interquartile range) of sources at z = 1.8 — 3.1. We find a 
statistically-significant difference between the measured redshift distributions in the 
two fields; the SXDF peaking at a slightly lower redshift (median z k, 2.2) than the 
LH (median z « 2.7), which we attribute to the noise-properties of the radio observa- 
tions. We demonstrate however that there could also be field-to-field variations that 
are consistent with the measured differences in the redshift distributions, and hence, 
that the incomplete area observed by SHADES with SCUBA, despite being the largest 
sub-mm survey to date, may still be too small to fully characterize the bright sub-mm 
galaxy population. Finally we present a brief comparison with the predicted, or as- 
sumed, redshift distributions of sub-mm galaxy formation and evolution models, and 
we derive the contribution of these SHADES sources and the general sub-mm galaxy 
population to the star formation-rate density at different epochs. 

(c) RAS 

Key words: surveys - galaxies: evolution - cosmology: miscellaneous - infrared: 
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1 INTRODUCTION 

The SCUBA HAlf Degree Survey (SHADES, Dunlop 2005, 
Mortier et al. 2005) was originally designed with the aim 
of characterizing the star-formation history (Hughes et al. 
2002) and clustering properties (van Kampen et al. 2005) 
of the bright-end of the luminous dust-enshrouded galaxy 
population. To achieve these goals we mapped two regions 
of the sky centered on the Lockman Hole East (LH) and 
Subaru XMM-Newton Deep Field (SXDF) with the Submil- 
limetre Common-User Bolometer Array (SCUBA, Holland 
et al. 1999). With a proposed la sensitivity of 2 mjy at 
850/im the complete survey was predicted to identify a sta- 
tistically robust sample of ~ 200 galaxies, with sufficient ra- 
dio to FIR ancillary data to help identify optical/IR counter- 
parts and derive spectroscopic/photometric redshifts. This 
redshift information is essential for determining the star for- 
mation and clustering properties for the whole population 
of ultraluminous dust-enshrouded galaxies. SCUBA was de- 
commissioned in mid-2005 having covered ~ 40 per cent of 
the originally-proposed area of the SHADES^. 

Paper I of this series (Mortier et al. 2005) describes the 
survey motivation, strategy and the philosophy adopted for 
the analysis. Paper II (Coppin et al. 2006) presents the cat- 
alogue and number counts derived from the 850/im sources. 
Paper III (Ivison et al. 2007) describes the identification of 
radio and mid-IR counterparts of these sources. This paper 
(IV) constructs the redshift distribution derived from the 
radio-mm-FIR photometry of the SHADES sources based 
on a compilation of the 850/im and 450/im SCUBA data 
(Coppin et al. 2006), 1.4GHz Very Large Array photometry 
(Ivison et al. 2007) and other previously published mm to 
FIR photometric observations towards these fields. A study 
of the mid-IR to optical properties of the SHADES popu- 
lation, and further constraints on the photometric redshifts 
of the sources, will be published elsewhere (Clements et al. 
2007, Dye et al. 2007, Serjeant et al. 2007 ). A spectroscopic 
study of a sub-sample of SHADES sources with identified 
optical/IR counterparts (Blain et al. 2007) will also provide 
an important comparison of spectroscopic and photometric 
redshifts. 

The cosmological parameters adopted throughout this 
paper are H = 71 kms" 1 Mpc" 1 , f2 M = 0.27, S1 A = 0.73. 



2 PHOTOMETRIC REDSHIFTS 

Despite having mapped only ~ 40 per cent of the planned 
0.5 sq. degree area, SHADES remains the largest extragalac- 
tic sub-mm survey to date. The difficulties of following-up 
such large areas at other wavelengths, and hence the inho- 
mogeneity of the multi-wavelength data, implies that the 
same photometric redshift technique cannot be applied to 
all sources. This section has been divided in two subsec- 
tions: the first (§ 2.1) deals with the consideration of 850/im 

1 The complete 1800 sq. arcmins SHADES area towards the LH 
and the SXDF has recently been surveyed at the JCMT at 1.1mm 
with AzTEC (Wilson et al. 2004), a continuum camera destined 
for the 50-m Large Millimetre Telescope (Serrano et al. 2006). 
These AzTEC data are currently being analysed and the results 
will be presented elsewhere. 



and 1.4GHz photometry which is available for all sources, 
and the use of the sub-mm-radio spectral-index as a diag- 
nostic of redshift; and the second (§ 2.2) describes the in- 
clusion of additional photometry at 70 to 450/tm which is 
sufficiently sensitive to place important constraints on the 
radio-mm-FIR photometric redshifts for only a few tens of 
sources. In both subsections we make a brief introduction 
to the techniques used, the estimated uncertainties found 
when comparing photometric and spectroscopic redshifts for 
similar sub-mm galaxies, the results from the application of 
the techniques to individual SHADES sources, and the com- 
bined redshift distributions derived for the entire SHADES 
population. 

2.1 1.4GHz/850/tm spectral index 

2.1.1 Techniques and accuracies 

One of the simplest redshift-indicators for the sub-mm 
galaxy population is that formed by the ratio of the flux 
densities at 1.4GHz and 850/im. These wavebands trace the 
tight correlation between radio continuum emission, which 
is dominated by synchrotron radiation from supernova rem- 
nants, and thermal emission from warm dust heated by 
young stars (Helou, Soifer & Rowan-Robinson 1985, Con- 
don 1992, Yun, Reddy & Condon 2001). This redshift in- 
dicator was systematically studied by Carilli & Yun (1999, 
2000), and has been subsequently revised for different sub- 
mm galaxy sub-populations (Dunne, Clements & Eales 2000, 
Rengarajan & Takeuchi 2001). The 1.4GHz to 850/im flux- 
density ratio, or a spectral index derived from it, increases 
monotonically with redshift, with some degeneracy due 
to the variety of radio synchrotron-slopes and mm dust- 
emissivity indices present in the ISM of those local galaxies 
used to define the relationship. Additionally there exists a 
level of degeneracy between the temperature of the dust gen- 
crating the rest-frame FIR luminosity (and hence sub-mm 
flux) and the redshift. Regardless, by adopting a library of 
local galaxy templates, and accepting the intrinsic disper- 
sion in their SEDs, the 1.4GHz to 850/im flux-density ratio 
still provides a crude but useful estimation of the redshift. 
This indicator becomes relatively insensitive to redshift be- 
yond z ~ 3, as the 850/im filter starts to sample the flatten- 
ing of the spectral energy distribution (SED) towards the 
rest-frame FIR peak, whilst still providing a powerful dis- 
criminant between low-redshift (z < 2) and high-redshift 
[z > 2) objects. 

We shall discuss the 1.4GHz/850/im spectral index fol- 
lowing two different prescriptions: (a) the single-template 
maximum-likelihood technique originally designed by Carilli 
& Yun (1999, 2000), denoted as Zp^ ot ; and (b) a maximum 
likelihood technique which simultaneously fits the 20 local 
templates of starbursts, ULIRGs and AGN used by Aretx- 
aga et al. (2003, 2005), denoted as Zp hot . 

The success of any photometric-redshift technique is 
measured by the accuracy with which it can predict the 
individual redshifts for a sample of representative galaxies 
with known redshifts, which have not been used to define the 
method. Aretxaga, Hughes & Dunlop (2006) have previously 
assessed the accuracy of the above two 1.4GHz/850/im pho- 
tometric redshift indicators. Based on this study, we show 
in figure 1 a comparison of spectroscopic and photometric- 
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Figure 1. (Top) Comparison of spectroscopic and photomet- 
ric redshifts derived from the 1.4GHz/850/im spectral index for a 
sample of 58 sub-mm galaxies with undisputed radio/optical/IR 
identifications, and spectroscopic redshifts derived from 2 or more 
lines (Aretxaga et al. 2006). The error bars represent 68% con- 
fidence intervals in the determination of the redshift. The r.m.s. 
scatter of the relation z spcc — z^ hot displayed is 0.8. (Bottom) 
Histogram distribution of the spectroscopic and photometric red- 
shifts represented in the top panel, which illustrates the success 
in recovering the redshift distribution of the sample. 



redshifts for 58 sub-mm/mm selected galaxies, comple- 
mented with a few objects selected at optical/FIR wave- 
lengths, which have published optical/IR or CO spectro- 
scopic redshifts and accompanying radio-FIR photometry. 
We will refer to this dataset as the 'comparison sample' 
hereafter. This comparison study shows that the Zp hot pre- 
scription has a mean accuracy Az =< \zp hot — z spcc \ >~ 0.65 
over the whole redshift interval, when one selects a robust 



sub-sample of objects with unambiguous optical/IR/radio 
counterparts and spectroscopic redshifts derived from the 
identification of two or more spectral lines. For the same ro- 
bust sample of objects, z^ ot has systematically larger errors, 
Az ~ 0.9. This sample does not include powerful radio-loud 
AGN, for which the template SEDs used in the photometric 
redshift analysis are not appropriate. The r.m.s. of the rela- 



tion is < (z phot - z s 



1/2, 



0.8. Restricting the analysis 



only to those galaxies with CO spectroscopic redshifts, the 
measured accuracy is Az « 0.6 and has an r.m.s of 0.8 at 
^ z 4. The precision degrades as the redshift increases, 
as expected from the 1.4GHz/850/im spectral index, which 
flattens beyond z — 3 (Carilli & Yun 2000), leading to a 
measured Az as 1.0 at 3 ^ z ^ 4. Using all objects with pub- 
lished photometry and spectroscopic redshifts, regardless of 
whether the associations that lead to the spectroscopic red- 
shift are unambiguous or not, the overall accuracy over the 
^ z ^ 4 regime degrades to Az ~ 0.8 (see Aretxaga et al. 
2006, figure 1). 

2.1.2 The redshifts of SHADES sources 

The radio counterparts adopted for the photometric redshift 
calculations of SHADES sources are the secure sample de- 
tected within 8 arcsecs of the sub-mm position (Ivison et 
al. 2007), with a chance-association probability between the 
radio and sub-mm source of P < 0.05. We have accepted 
some additional counterparts when a robustly-detected ra- 
dio source is still within 10 arcsec of the sub-mm centroid 
and, additionally, a 24/mi counterpart is associated with 
this radio identification. These extra radio counterparts are 
marked in the notes provided for each sub-mm source (see 
tables [1] and [2j| , where we have calculated the correspond- 
ing P-value of the radio association, which remains lower 
than 0.08. The 34 radio sources adopted as counterparts of 
SHADES galaxies in the LH field, and the 35 radio sources 
in the SXDF have a combined chance association P ~ 1.6, 
and thus we expect to have incorrectly associated ~ 1 of the 
SHADES sub-mm sources with a projected radio source. 

For both techniques the error bars of the photometric 
redshifts were derived by bootstrapping on the reported pho- 
tometric and calibration errors (Coppin et al. 2006, Ivison 
et al. 2007), and are defined as the 68% confidence interval 
of the resulting redshift probability distribution. The pho- 
tometric error distributions used for the 850/im photometry 
were derived by de-boosting the measured flux densities of 
the SCUBA sources. A de-boosting correction is necessary 
to provide a more accurate estimate of the flux of low S/N 
blank-field sources in sub-mm surveys, where the counts are 
typically very steep and faint galaxies can be statistically 
boosted above the nominal detection threshold. These errors 
are often non-Gaussian (see figure 5 in Coppin et al 2006). 
The 1.4GHz flux densities do not need to be de-boosted, 
since this correction is dependent on the area defined by the 
search radius in identifying the source. In the case of finding 
associations within <J 8 arcsec radius around a known ob- 
ject, this is negligible. The error distributions for the 1.4GHz 
flux densities were assumed to be Gaussian. In the case of 
z photi the error estimated by Carilli & Yun (2000), to allow 
for a difference in templates, is added in quadrature to the 
errors derived by bootstrapping the photometry. 

The probability distribution calculated for each source 
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Figure 2. (Top) Spectral energy distribution (SED) of 
SXDF850.8, where the black squares mark the detections at 
850/xm and 1.4GHz used in the photometric redshift calculation. 
Error bars are la, and the arrow at 450/mi marks the 3cr up- 
per limit derived from our maps. For reference, the SED tem- 
plates used in the photometric redshift calculation are shifted 
to Zp hot = 2.6 and scaled to maximize the likelihood function 
of detections and upper limit through survival analysis (Isobe, 
Feigelson & Nelson 1986), and are represented as lines. All the 
SEDs are compatible within the 3cr error-bars of the photometry 
of the source. (Bottom) Probability distribution for SXDF850.8 
derived for the z^ hot solution, using only the 850//m and 1.4GHz 
photometry. 

in this manner typically have a single peak, which broadens 
as the most probable redshift of the source increases. Fig- 
ure 2 shows an example of a typical solution derived from 
the use of 850/xm and 1.4GHz photometry. 

2.1.3 Redshift distribution of the SHADES population 

We have assembled the best estimates of photometric red- 
shift for each source in the LH and SXDF SHADES fields 
(figure 3) . These are identified by the modes of the individ- 



Figure 3. Histogram of modes of the photometric redshift dis- 
tributions of SHADES galaxies derived from the 1.4GHz/850/^m 
spectral index. The (black) thick solid-line (shown in the upper 
and lower-panels) represents the distribution of modes for the 69 
galaxies that have been detected at both 850/im and 1.4GHz. In 
the upper-panel the (blue) thin dotted-line and (red) thin dash- 
dotted line represent the redshift distributions in the LH and 
SXDF fields respectively. In the lower-panel the black dashed-line 
(a) and grey dashed-line (b) show the redshift distributions for 
the full SHADES catalogue, including the 51 sub- mm galaxies un- 
detected at 1.4GHz. Those SCUBA galaxies with non-detections 
in the radio are distributed in one of two ways that bracket the 
range of reasonable options: (a) with equal probability between 
their calculated lower 90% confidence limits and z = 5; and al- 
ternatively, (b) between their lower limits and z = 2, or only at 
their lower redshift-limits in the cases that these lie at z > 2. 

ual probability distributions, Zp hot , which, by definition, are 
the redshifts with the highest probability values. 

The galaxies that are not detected with confidence at 
radio- wavelengths, i.e. 26 out of the 60 galaxies in the LH, 
and 25 out of the 60 in the SXDF, have very flat individual 
redshift probability distributions in our computations, and 
hence we quote only the lower-limits to their redshift, which 
are defined as their 90% lower confidence-limits (tables 1 
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and 2). We incorporate these objects into the population 
distribution by adding flat probability distributions between 
their calculated lower 90% confidence limits and z = 5, and 
alternatively between their calculated lower 90% confidence 
limits and z — 2, or only at their lower redshift-limits if these 
indicate z > 2. These two alternative priors illustrate how 
the resulting redshift distributions (that include SHADES 
galaxies without radio detections) can be biased high and 
low. 

Figure shows the final photometric redshift distribu- 
tion using the 1.4GHz/850/xm spectral index, both for the 
full SHADES sample and for the LH and SXDF fields sep- 
arately. The redshift distribution of SXDF sources peaks 
at slightly lower redshifts (median z ft! 2.2) compared to 
the distribution of LH sources (median z w 2.6). The low- 
redshift tail (z < 1.5) is also slightly more prominent in the 
SXDF than in the LH (6 vs. 1 sources among the radio- 
detected sample). The difference in shape of the two dis- 
tributions of radio-detected sources can be measured using 
the two-tailed Kolmogorov-Smirnov (K-S) test, which gives 
a 3 per-cent chance that they are drawn from the same par- 
ent distribution. Furthermore, the mean-redshifts of the two 
distributions are significantly different at the 99.7 per cent 
level, according to a Mann- Whitney [/-test. This difference 
can be attributed to differences in the noise levels of the 
radio maps, although some intrinsic variations between the 
fields are expected (see § 3). 

For those objects with more than one redshift estimate 
(due to the ambiguity in their counterparts), we have pro- 
duced alternative population distributions, with or without 
their inclusion, and with different combinations of possible 
counterparts. The results do not significantly change the fi- 
nal combined distribution. All figures presented in this paper 
include the primary radio association if there are multiple 
options. 

2.2 Radio-mm-FIR SED analysis 

The SHADES fields have been targeted by other FIR/sub- 
mm/mm and radio surveys. In this subsection we describe 
the extra constraints on the photometric redshifts that can 
be derived from these additional complementary data for a 
few tens of SHADES sources. 



2.2.1 Techniques and accuracies 

Photometric redshifts with modest precisions (Az ~ 0.3 
to 0.5) have been obtained in the past few years using a 
combination of spectral indices between the radio and mm- 
wavelength regimes and the FIR spectral peak. This infor- 
mation has been exploited by several groups using a wide ar- 
ray of fitting-techniques and SEDs (e.g. Yun & Carilli 2002, 
Hughes et al. 2002, Aretxaga et al. 2003, 2005, Wiklind 2003, 
Hunt & Maiolino 2005, Laurent et al. 2006). There remain, 
however, degeneracies imposed by the choice of multiple 
SED templates, with FIR emission peaks distributed over 
a range of wavelengths, which can limit the precision of the 
derived redshifts (e.g. Blain, Barnard & Chapman 2003). 
We have previously developed a radio-mm-FIR technique 
based on Monte-Carlo simulations, that take into account 
constraining prior information such as the number counts 



of sub-mm galaxies, the favoured luminosity /density evolu- 
tion up to z ~ 2, and the lensing amplification of a certain 
field (Hughes et al. 2002, Aretxaga et al. 2003, 2005). We 
only offer a brief summary of this technique here. A cata- 
logue of 60/im luminosities and redshifts for mock galaxies 
is generated from an evolutionary model for the 60^tm lumi- 
nosity function that fits the observed 850/im number-counts 
(e.g. luminosity evolution oc (1 + z) 3 for z <, 2, and no 
evolution at z > 2) and covers a simulated area of 10 sq. 
deg. Template SEDs are drawn at random, without regard 
to their intrinsic luminosity, from a library of 20 local star- 
bursts, ULIRGs and AGN, to provide FIR-radio colours for 
the mock galaxies. The SEDs cover a wide-range of FIR 
luminosities (9.0 < logZ/FiR,/ 'Lq < 12.3) and temperatures 
(25 < T/K < 65). The flux densities of the mock galaxies 
include both photometric and calibration errors, consistent 
with the quality of the observational data for each sub-mm 
galaxy detected in a particular survey. We reject from the 
catalogue those mock galaxies that do not respect the detec- 
tion thresholds and upper-limits of the particular sub-mm 
galaxy under analysis. The redshift probability distribution 
of an individual sub-mm galaxy is then calculated as the nor- 
malized distribution of the redshifts of the mock galaxies in 
the reduced catalogue, weighted by the likelihood of identi- 
fying the colours and flux densities of each mock galaxy with 
those of the sub-mm galaxy in question. This technique will 
be denoted z^o t in the discussion that follows. 

The validity of the results derived from this technique 
is limited by the assumption that the SEDs of high-z sub- 
mm galaxies are similar to the local analogues, adopted as 
templates, which are scaled in luminosity and shifted in 
redshift. While this might seem a naive approach, all the 
templates used in the calculations that follow offer a good 
description of the radio-mm-FIR photometry of SCUBA 
galaxies, including 350/im observations (Laurent et al. 2006, 
Kovacs et al 2006), with known spectroscopic redshifts and 
un-ambiguous multi-wavelength counterparts (Aretxaga et 
al. 2005). There are however a few sub-mm galaxies which do 
not match any of the templates we use in this paper at their 
published redshifts (see figure 4 in Aretxaga et al. 2006). 
In these examples their radio emission is higher than that 
implied by the radio-FIR correlation, possibly due to ac- 
cretion activity, or their FIR emission peaks at wavelengths 
longer than those of the templates used in this study at 
the adopted redshift. We describe the redshift solutions for 
these galaxies as 'catastrophic' and this might be indica- 
tive of incompleteness in the library of SED templates used 
as analogues. There is still sufficient debate in the litera- 
ture, however, about the nature and/or the ambiguity of 
the multi-wavelength counterparts to these sub-mm sources, 
from which the redshifts are derived, to justify their exclu- 
sion from a robust comparison sample (see Aretxaga et al. 
2005, 2006, Laurent et al. 2006, Kovacs et al. 2006 for a 
detailed complementary discussion on these galaxies). 

In order to estimate the accuracy of the z^, t technique 
we use the full SED information of a robust sub-sample of 
11 sub-mm galaxies, out of the comparison sample of 58 
galaxies considered in §2.1.1, which have detections in three 
or more bands. Furthermore, the same galaxies have undis- 
puted identifications of their optical/IR/radio counterparts 
and spectroscopic redshifts derived from the measurement 
of two or more spectral lines. We derive a mean accuracy 
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Figure 4. (Top) Comparison of spectroscopic and photomet- 
ric redshifts derived from the full radio-FIR SED for a sample of 
11 sub- mm galaxies with at least 3 robust detections at different 
wavelengths. This sample has undisputed radio/optical/IR coun- 
terparts associated with the sub-mm galaxies, and spectroscopic 
redshifts derived from 2 or more lines. The relationship has an 
r.m.s. of 0.25. (Bottom) Comparison of the distributions of the 
spectroscopic and photometric redshifts represented in the top 
panel. 



for this sub-sample of Az 

,MC 



r.m.s. < (Zphot — Zspcc) 2 > 1 ^ 2 « 0.25 over the whole redshift 
interval (figure 4). Using all objects with published photom- 
etry, regardless of whether the spectroscopic redshift derived 
from the optical associations is ambiguous or not, the overall 
accuracy over the ^ z ^ 4 regime degrades to Az w 0.55, 
with an r.m.s. of 0.80 (see figure 3 in Aretxaga et al. 2006). 
A few significant outliers which remain in the correlation 
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are discussed by Aretxaga et al. (2005) and Kovacs et al. 
(2006). Within the small sub-sample of study, the accuracy 
is independent of redshift. 

If we restrict the use of photometry to 450/xm upper 
limits combined with the 1.4GHz and 850/im detections 
for the comparison sample of galaxies (adopting simulated 
450^tm upper limits, when necessary, to mimic a shallower 
survey at this wavelength), we find a mean accuracy of 
Az ~ 0.55 and an r.m.s. of 0.7 for the robust sample. This 
result is especially relevant for the photometric redshift cal- 
culations of the SHADES sources in £12.2.21 the majority of 
which have similarly sparsely-sampled photometry. Consid- 
ering only the complete sample with robust and tentative 
spectroscopic redshifts the mean accuracy degrades slightly 
to Az » 0.65, and the r.m.s. scatter is 0.90. 



2.2.2 The redshifts of SHADES sources 

Table 3 summarizes the most recent photometric redshifts 
calculated with the Monte Carlo technique for SHADES 
sources with additional photometry published in the litera- 
ture. In contrast, and for completeness, tables 4 and 5 list 
the photometric redshifts derived only from the combina- 
tion of the SHADES 450/850/mi and the 1.4GHz photome- 
try using two approaches: the Monte Carlo based technique, 
z phot described above, and a non-prior maximum likelihood 
fit to the same 20 SEDs used for the first method that in- 
cludes survival analysis (Isobe, Feigelson & Nelson 1986) to 
incorporate the non-detections into the maximum likelihood 
formalism, Zp* ot . This second technique is introduced to pro- 
vide a comparison of how the priors affect the redshift esti- 
mation of the sources and the final combined redshift distri- 
bution of SHADES galaxies. While most of the photometric 
redshifts derived from the two methods are similar, the pure 
survival analysis produces a few high-z catastrophic results 
in the robust comparison sample (2 out of 11). The overall 
reliability of the maximum likelihood technique is Az ~ 0.7. 
These high-z catastrophic solutions get suppressed by the 
MC technique due to the weighting priors that disfavour 
high-z solutions for these sources, since, if they were typi- 
cal of the sub-mm population, they would overproduce the 
850/im number counts under the assumed luminosity evo- 
lution model. Although we give the values of photometric 
redshifts with and without priors in tables 4 and 5, we will 
now continue the analysis of the complete SHADES sample 
using only the MC solutions, since they have been shown to 
perform better against the comparison sample. 

For 5 sources in table 5, SXDF850.5, 21, 28, 77 and 
119, we also include complementary photometry at 70 and 
160^im from the Spitzer Legacy Survey SWIRE (Lonsdale et 
al. 2003, Surance et al. 2007) that are used to derive mid-IR 
counterparts to the SHADES sources (Clements et al. 2007). 
The remainder of the SHADES sources are not significantly 
detected (> 4a) in the Spitzer catalogues, and the noise 
properties of the SWIRE maps, providing 3<r upper-limits 
of ~23 mjy and 160 mjy at 70 and 160/xm, respectively 
(Afonso-Luis et al. 2007), do not further constrain the pho- 
tometric redshifts. 

Figure 5 shows the single-peaked redshift probability 
distribution derived for one of the sources that have the 
most complete photometric data. Although the majority of 
the sources show similar probability distributions, there are 
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Figure 5. (Top) Spectral energy distribution (SED) of LH850.1, 
where the black squares mark detections, with la error bars, and 
arrows indicate 3cr upper limits. For reference, the SED tem- 
plates used in the photometric redshift calculation are shifted 
to z^ c = 2.4 and scaled to maximize the likelihood function 
of detections and upper limit through survival analysis (Isobe 
et al. 1986), and are represented as lines. The SEDs compatible 
within the 3<r error- bars of the photometry of the source are repre- 
sented in darker (blue) lines. (Bottom) Probability distribution 
of LH850.1 derived for the z^ c solution. 

sometimes secondary peaks (see examples in Aretxaga et al. 
2003). Nevertheless, it is always the primary redshift-peak 
that defines the solutions given in tables 3,4,5. 

2.2.3 Redshift distribution of the SHADES population 

Figure 6 shows our final photometrically-derived redshift 
distribution for the SHADES sources, using our best avail- 
able estimate for the redshift of each source (i.e. taken 
from table 3 for those sources with the most complete pho- 
tometry, and from tables 4 and 5 for the remainder) 

The distribution of radio-identified SHADES sources 
clearly peaks in the bin z m 2.0 — 2.5, with a 50 per cent 



interquartile interval z ~ 1.8 — 3.1. We incorporate the non- 
radio detected sources (lower panel in figure 6) into the pop- 
ulation distribution in two alternative ways, to serve as ex- 
amples of how much these sources could alter the final popu- 
lation distribution: (a) approximating their individual prob- 
ability distributions as flat distributions between their lower 
90% confidence limit and z — 5; and (b) as flat distributions 
between their lower 90% confidence limit and z — 2, or at 
their lower limit if this lies at z > 2. Solution (a) is actually 
derived from the adopted non-informative (flat) prior for the 
photometric-redshift calculations. This creates a high-z tail 
which is a reflection of the adopted range for the flat redshift 
distributions which are unconstrained by the photometry. 
Solution (b) is biased against high-z, by imposing a maxi- 
mum redshift for the radio-undetected sample which is lower 
than the redshift of the peak of the radio-detected sample. 
This radio-undetected sample could be composed of colder 
sub- mm galaxies than those found in the template library, or 
they could have the same template shapes as those adopted 
in the photometric redshift analysis, and still be undetected 
at the depth of the present radio surveys. Regardless, in 
these alternative solutions, the mode of the population re- 
mains at z « 2.0 — 2.5, with at least 50 per cent of the 
galaxies in the interquartile range 1.6 ^ z ^ 3.4. 

In order to consider the effect of the objects with more 
than one redshift estimate (due to ambiguity in their radio- 
counterparts), we have produced alternative population dis- 
tributions. For instance, figure 6 shows the combination of 
the first entries for each source in tables 3, 4, 5. The in- 
troduction of the second tabulated values instead of the 
first ones, for those sources with ambiguous associated pho- 
tometry, produces an alternative distribution which is in- 
distinguishable (with a 99.96 per-cent probability), via a 
Kolmogorov-Smirnov (K-S) test, from the one represented 
here. 

As in the case of the photometric redshift distribution 
derived from the 1.4GHz/850/im spectral index, the dis- 
tribution of redshifts derived from the full SED analysis 
of SXDF sources peaks at slightly lower redshifts (median 
z ~ 2.2) than that of LH sources (median z « 2.7), and its 
low-redshift tail (z < 1.5) is also more prominent. These dif- 
ferences in the distributions are statistically-significant, as 
indicated by a K-S test at a level of 98.9%. A Mann- Whitney 
U-test shows that their mean-redshifts differ at the 99.997% 
level. 



3 DISCUSSION 

3.1 Redshift distribution 

SHADES was designed with the objective of constraining 
the redshift distribution and clustering properties of the sub- 
mm galaxy population, an exercise which van Kampen et al. 
(2005) demonstrated could discriminate between galaxy for- 
mation models. With ~40 per cent of the survey completed 
before SCUBA was de-commissioned in the summer 2005, 
SHADES has provided 120 robust sources. The radio-mm- 
FIR photometry assembled for the survey favours a nearly 
Gaussian redshift distribution of the population peaking at 
z ~ 2.0—2.5, albeit still with the possibility of a high-redshift 
tail remaining. 
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Figure 6. Histogram of the modes of the photometric redshifts of SHADES galaxies based on all available radio-mm-FIR photometry 
(provided in table 3, or otherwise in tables 4 and 5). Upper Panel: The thick solid-line (also shown in the lower-panel) represents 
the z^f^ distribution of modes for those 70 SCUBA galaxies that have been detected in at least two bands, of which the thin (blue) 
dotted-line and thin (red) dash-dotted line represent the distributions defined by the LH and SXDF sources, respectively. For comparison, 
we plot the Zpj^ ot solutions for the full distribution of LH and SXDF (dash-3-dot line). Lower Panel: The black dashed-line (a) and 
grey dashed-line (b) show the redshift distributions for the 70 SCUBA galaxies detected in at least two bands, plus an additional 50 
SCUBA galaxies detected only at 850/^m. These latter sources are distributed in one of two ways that bracket the range of reasonable 
options: (a) with equal probability between their calculated lower 90% confidence limits and z = 5; and alternatively, (b) between their 
lower limits and z = 2, or only at their lower redshift-limits in the cases that these lie at z > 2. 



The photometric redshift distribution of the radio- 
detected sub-mm galaxies is qualitatively similar to the op- 
tical spectroscopic redshift distribution published by Chap- 
man et al. (2003, 2005) who followed-up a sample of sub- 
mm galaxies derived from various surveys. The agreement 



is perhaps not surprising, given that the photometric red- 
shifts of the comparison sample have shown a relatively good 
agreement with the spectroscopic redshifts published in the 
literature (Aretxaga et al. 2006, §2.1.1, 2.2.1). Furthermore, 
the majority of the sub-mm sources that have spectroscopic 
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redshifts are drawn from SCUBA surveys of similar depths 
to SHADES. 

The high-redshift correction applied to the measured 
spectroscopic redshift distribution, suggested by Chapman 
et al. (2005) to account for the bias introduced by non- 
detection of the higher-redshift radio counterparts that pro- 
vide candidates for optical spectroscopic follow-up, also falls 
within the range of photometric redshift estimations we have 
derived for SHADES sources that are not detected at radio- 
wavelengths. These sources provide the high-redshift (z > 3) 
tail of figure 6, and could in fact be placed anywhere above 
z ~ 1.0, even producing secondary peaks. 

Our calculations do not support the existence of a sub- 
stantial low-redshift (z < 1.5) tail within the luminous sub- 
mm radio-detected population sampled by SHADES. At first 
sight this might appear to be in conflict with the results of 
Pope et al. (2005, 2006) who found that ~ 30% of the sub- 
mm sources found in the SCUBA imaging of the GOODS- 
North field may lie at z < 1.5. If the GOODS-North 850/wn 
catalogue is restricted to sources with de-boosted flux densi- 
ties 5850 > 3mJy, however, then the proportion of robustly 
identified sub- mm galaxies which lie at 2 < 1.5 drops to 
6%. This is entirely consistent with the results found here 
for SHADES galaxies. Thus, these results may be provid- 
ing further evidence that the peak of the redshift distribu- 
tion of sub-mm sources is positively correlated with sub- 
mm flux-density/luminosity, consistent with the apparently 
anti-hierarchical nature of star-formation history reported in 
several other recent studies (e.g. Heavens et al. 2004). Fur- 
thermore field-to-field variations in the spatial distribution 
of the large-scale structure can provide a simple explana- 
tion for the differences between the redshift distributions 
of sub-mm sources derived from the individually-mapped 
contiguous-areas (typically < 0.2 sq. degrees) taken from 
the current generation of SCUBA surveys. 

Our analysis also suggests that only a modest fraction 
of sub-mm galaxies could be hiding in the optical redshift- 
desert atz^s 1.5 — 1.8 during spectroscopic searches for 
SHADES sources with robust radio counterparts. The pho- 
tometric redshift probability density distributions of radio- 
detected SHADES sources using the 1.4GHz/850pm index 
or the full radio-mm-FIR SED information contain ~ 15 % 
and ~ 10% of sources in this redshift desert regime, respec- 
tively. 

The difference between the redshift distribution of sub- 
mm sources in the SXDF and LH fields is entirely consis- 
tent with the different properties of the 1.4-GHz maps as 
discussed by Ivison et al. (2007). The LH data have a higher 
VLA resolution than those in SXDF, and the LH data are 
also deeper, although the coverage is less uniform. There 
is clearly potential for systematic differences between radio 
measurements in the LH and SXDF. For an extended source 
in the LH (of which there are several - Ivison et al. 2002), a 
larger fraction of emission on scales larger than the synthe- 
sized beam will be resolved away than for similar cases in 
SXDF. Moreover, the LH data will suffer greater significant 
bandwidth smearing and, although the appropriate correc- 
tion has been made to the measured flux densities, some 
faint sources will be lost below the radio-detection threshold 
and may receive misleadingly low flux-density limits. These 
effects can be viewed as a systematic flux calibration offset 
with consequences as severe as those encountered in opti- 



cal/infrared photometric-redshift estimation. While random 
1.4GHz calibration uncertainties of 5 per cent have been ac- 
counted for in the estimation of the photometric redshifts, a 
systematic flux-density offset could shift the redshift distri- 
bution significantly. In order to explore this possibility, we 
have applied a 10 per cent flux increase to the LH photom- 
etry and recalculated the photometric redshifts. The com- 
bined redshift distribution shifts its peak by ~ —0.25, and 
consequently the mean values of the SXDF and LH distri- 
butions are more consistent, increasing from 0.3 per cent 
(§2.1.3) to a 7 per-cent probability, according to a Mann- 
Whitney [/-test. Some intrinsic variation on the distribution 
of redshifts between the fields is however to be expected (see 
below). 

3.1.1 Comparison of the SHADES redshift distribution 
with galaxy formation models 

Van Kampen et al. (2005) studied four different galaxy for- 
mation models that yielded different redshift distributions 
and clustering properties for the sub-mm population ex- 
pected to be found in a survey of the depth and area cov- 
ered by SHADES: (a) a hydrodynamical model (Muanwong 
et al. 2002), that follows the evolution of dark-matter, gas, 
star-like particles and galaxy fragments, that has been cou- 
pled with the analytical form for redshift distribution of 
Baugh, Cole & Frenk (1996); (/3) a simple merger model 
that identifies sub-mm galaxies with major mergers of mas- 
sive galaxies; (7) a phenomenological model (van Kampen 

2004) , which is based on TV-body simulations that identify 
the sites of major-mergers and has two modes of star for- 
mation, quiescent and bursting; and (S) a stable clustering 
model (Gaztanaga & Hughes 2001). Figure 7 represents the 
theoretical redshift distributions of SHADES galaxies found 
in these models. This figure has been complemented with (e) 
a semi-analytical model for the joint formation and evolu- 
tion of spheroids and QSOs (Granato et al. 2004, Silva et al. 

2005) ; and (£) an alternative semi-analytic model of galaxy 
formation for sub-mm galaxies (Baugh et al. 2005). 

Furthermore, to enable a more accurate discrimination 
between the above predictions, all the galaxy formation 
models in figure 7 account for the incompleteness of sources 
in the SHADES catalogue (Coppin et al. 2006). The models 
have also been convolved with a representative radio-mm- 
FIR photometric precision of a ~ 0.4, which is intermedi- 
ate between the measured uncertainties derived for the two 
techniques used in this paper. 

We have made a comparison, via a K-S test, of the ob- 
served redshift probability density distributions with those 
predicted from the above models. In each case a K-S statistic 
has been calculated that accommodates the la uncertainty 
in the median redshift of the models due to field-to-field vari- 
ations. A study of 25 simulations made for each of the four 
models analyzed by van Kampen et al. (2005) demonstrate 
that the mean redshift of ~ 60 SHADES-like sub-mm galax- 
ies varies by Sz (r.m.s.) « 0.25 — 0.55. In part these shifts 
can be explained by Poisson noise (estimated a ~ 0.1 — 0.2 
from the simulations). The models show, however, that there 
could also be a significant component in the field-to-field 
variations that arises from intrinsic redshift differences due 
to varying amounts of groups or proto-clusters of galaxies 
along the line-of sight. Thus the differences found between 



© RAS, MNRAS 000, 



10 Aretxaga et al. 




MC 

phot 



Figure 7. Probability density of the combined redshift distribution of SHADES galaxies (thin black solid-line, and thin dashed-lines ((a) 
and (b)), as described in figure 6). These are compared with the redshift distributions of six galaxy formation models (thick coloured- 
lines), degraded with a a z = 0.4 to provide representative redshift-uncertainties: (a) the hydrodynamical model of Muanwong et al. 
(2002) coupled with the analytical form for redshift distribution of Baugh, Cole & Frenk (1996); (/3) the simple merger model of van 
Kampen et al. (2005); (7) the phenomenological model of van Kampen (2004) and van Kampen et al. (2005); (5) the stable clustering 
model of Gaztanaga & Hughes (2001); (e) the semi-analytical model for the joint formation and evolution of spheroids and QSOs of Silva 
et al. (2005); and (f) the semi-analytic model for galaxy formation of Baugh et al. (2005). 



the LH and SXDF areas, and between these and smaller, 
deeper surveys like GOODS-N, could be partially explained 
by this effect. 

The results of the K-S test suggest that only model (e) is 
close to being formally acceptable, with an 87% probability 
for the model to agree with the measured probability density 
distribution that includes SHADES sources with and with- 
out radio-detections according to solution (a). With only a 
small shift (Sz ~ —0.3) in the distribution, model (e) also 
qualitatively reproduces (~ 60% probability of similarity) 
the photometric-redshift distribution of the radio-detected 
SHADES galaxies. 

The SHADES sources in our analysis that are not de- 
tected at radio wavelengths have very flat redshift probabil- 
ity distributions, which simply places them at z ^ 1.0, and 
hence these SHADES sources could also produce a secondary 
peak in the redshift distribution. In the ranking of similar- 
ities of measured and model-distributions, models (a) and 
(7), ~ 45% probability, have double peaks and are broader 
than the observed distributions. A different prior, that opti- 
mizes the redshift-distribution of the SHADES sources with- 
out radio-detections, could bring them closer to a level of 
formal-acceptance. Finally models (S), (£) and (/3) are all 
rejected with probabilities of < 2% of being consistent with 
the range of solutions depicted in figure 7. 



3.2 The FIR luminosity of SHADES sources 

The catalogues of redshifts presented in tables 1 to 5 are 
an initial step towards characterizing the FIR luminosities 
and star formation rates of the SHADES population. The 
available photometry in the FIR peak regime (70-450/^m), 
however, is not deep enough to fully constrain the SEDs 
of most SHADES sources at these wavelengths. One viable 
approach is to use the 20 SEDs in our local template cat- 
alogue to derive the corresponding FIR luminosities from 
the 850/im flux densities, bearing in mind that the lack of 
constraints at short wavelengths will dominate the errors in 
luminosity estimation over those of redshift (e.g. Hughes et 
al. 2002). Alternatively, one could use the 1.4GHz radio flux 
density to deduce FIR luminosities via the radio-FIR lu- 
minosity correlation that characterizes the sub-mm galaxy 
population, since this now has been extended to 2 ~ 0.5 — 4 
(Kovacs et al. 2006). This latter approach has the advan- 
tage of providing mean FIR luminosities which are accurate 
for the bulk of the population, reducing the uncertainties 
in luminosity primarily to the accuracy of the photometric 
redshifts. However, the normalization of the relation might 
be shifted from the local IRAS correlation, and this could 
affect the FIR luminosities derived, and the comparison of 
these to nearby galaxies. 

Regardless of this complication, the observed 1.4GHz 
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Figure 8. Distribution of FIR luminosities of the 120 SCUBA 
galaxies in the SHADES catalogue. The solid-line represents the 
total distribution of radio-detected sub-mm sources in the two 
fields towards the LH (thin blue dotted-line) and SXDF (thin red 
dashed-dotted line) based on the FIR-radio correlation of sub-mm 
galaxies (Kovacs et al. 2006). The dashed line shows the complete 
distribution of FIR luminosities for all SHADES sources, includ- 
ing the 51 non-radio-detected sub-mm sources, where the FIR 
luminosities of the latter sources have been derived from a single 
T = 35K, /3 = 1.5 grey-body scaled to the observed 850/im flux 
density, and their redshifts have been selected at random between 
their lower 90%-confidence redshift-limits and, arbitrarily, 2 = 5. 



flux densities have been converted to rest-frame 1.4GHz flux 
densities using a mean synchrotron radio slope of index a = 
—0.7, and the monochromatic 1.4GHz luminosity has been 
inferred using the photometric redshift solution for each 
source. This is converted to FIR luminosity using the lin- 
ear relationship log(L F m/Li.4GHa/4.52THz) = 2.14 ± 0.07 
(Kovacs et al. 2006). For each source we have considered 
the effect of the uncertainties in redshift, 1.4GHz flux and 
the reported scatter in the sub-mm galaxy FIR-radio cor- 
relation (Kovacs et al. 2006) by bootstrapping 1000 times 
on the measured errors. For the 69 radio-detected SHADES 
sources in the LH and SXDF fields the median FIR luminos- 
ity is 2.6 x 1O 12 L0, with a high luminosity tail that extends 
to 1 x 1O 13 L0 (see figure|8}. The distribution of luminosities 
for sources in both the LH and SXDF fields are similar. 

The effect of non-radio-detected sources, which have 
very unconstrained and possibly high redshifts (§2.1), has 
also been considered by scaling a simple grey-body of tem- 
perature T — 35 K and emissivity index /3 = 1.5, the aver- 
age of a parametrized SED of the short sub-mm wavelength 
detected SCUBA galaxies that define the radio-FIR corre- 
lation at z ~ 1 — 3 (Kovacs et al. 2006), to the observed 
850/im flux density. The errors in the 850/im flux density 
have been taken into account by bootstrapping on the in- 
ferred de-boosted distributions of 850/im flux densities for 
each source, and the redshift has been selected at random 
between their 90 per cent lower-limits (tables 4, 5) and, ar- 
bitrarily, z — 5. The resulting combined distribution has a 
median FIR luminosity of 2.6 x 1O 12 L0, the same as the 
distribution of radio-detected sources. 

If one returns to the alternative approach of using 



the combination of 20 SEDs with the de-boosted 850/im 
flux densities to derive the FIR luminosities, one derives a 
broader luminosity distribution than that depicted in fig- 
ure 8, reflecting the wide variety of acceptable SED tem- 
plates, with a median luminosity that is increased by ~ 
40 per cent. 

3.3 The star formation-rate history derived from 
SHADES galaxies 

The evolution of the global star formation rate (SFR) den- 
sity traced by SHADES sources is shown in figure [9] The 
conversion from FIR luminosity to star formation rate was 
performed using a constant of 5 x 10 9 L0/(M0yr _1 ), which 
is constrained to ~ ±30 per cent uncertainty (Kennicutt 
1998). We have multiplied the contribution of each source 
to the star formation density by the inverse of the SHADES 
survey completeness function at the appropriate flux den- 
sity (Coppin et al. 2006). Redshift-space was binned into six 
intervals, and a Monte Carlo was performed to assign each 
galaxy to a redshift bin according to its expected photomet- 
ric redshift error. The FIR luminosities have been estimated 
from the radio-FIR correlation or with a single SED which is 
considered to be representative of the sub-mm galaxy pop- 
ulation (as in §3.2). The error bars in SFR density are the 
result of the uncertainties in photometry, SEDs and redshift, 
and are computed as the standard deviation traced by 1000 
Monte Carlo simulations. 

The evolution in the SFR density traced by the radio- 
detected sources shows a clear peak at z ~ 2.5 and a slow 
decline at both low and high redshifts. A good descrip- 
tion of the SFR density traced by this population is given 
by psf « 0.35 exp [-0.5(z - 2.8) 2 /(0.8 2 )] . The contribu- 
tion of non-radio-detected SHADES sources to the SFR his- 
tory may be significant at high-redshifts compared to the 
radio-detected SHADES sources, as indicated in figure 9 by 
the empty black diamonds (solution (a) in figure 6). Since 
the redshifts of these sources are not well constrained, this 
should be considered as only a possible evolutionary his- 
tory, awaiting confirmation by better multi- wavelength data 
to improve the constraints on the redshifts of the radio- 
undetected SHADES sources. The effect of placing all the 
radio-undetected sources between their lower redshift-limits 
and z = 2 or at their lower limits if they are at z > 2 (solu- 
tion (b) in figure 6) is shown by the grey empty diamonds 
in figure 9, which understandably broadens the peak of star 
formation density to lower- z. 

The levels of star formation deduced for SHADES 
sources are consistent with those derived for other samples of 
radio-detected sub-mm galaxies at z ^ 3.5 (Chapman et al. 
2005). At higher-redshifts, however, the SFR densities im- 
plied in this photometric-redshift study exceed the extrap- 
olations of Chapman et al. Eventually spectroscopic mea- 
surements of SHADES sources at z > 3.5, via millimetre- 
wavelength observations of molecular CO-lines, or optical 
spectroscopy with increased sensitivity, will provide a defini- 
tive measurement of the SFR density of obscured galaxies 
in the high-redshift Universe. 

In order to estimate the contribution to the global star 
formation rate density of 850/im sources that are fainter 
than those detected by SHADES, we have adopted the 60/im 
luminosity function with a pure luminosity evolution that 



© RAS, MNRAS 000, 



12 Aretxaga et al. 



0.100 



u 

Q_ 



0.010 



0.001 








2 



3 
z 



4 



6 



Figure 9. Evolution of the global star formation rate density (SFRD) for different samples of galaxies. We show the SFRD for radio- 
detected SHADES sources (solid diamonds) and for all SHADES galaxies, considering both the radio-detected and radio-undetected 
sources (open diamonds). The redshift probability distributions for the radio-undetected sources have been assumed to be flat between 
their lower redshift-limits (tables 4, 5) and 2 = 2 (grey open diamonds) and z = 5 (black open diamonds). The redshift error-bars 
indicate the width of the redshift bins. The error-bars in p$F are a combination of the uncertainty in the photometry, the recovery of the 
luminosity of each source due to SED differences, and the uncertainty in redshift that divides the contribution of sources across several 
redshift bins (see §3.3). The empty-square grey/black symbols represent the SFRD traced by ultraluminous 850/mi-selected starbursts, 
by correcting the SFRDs derived for SHADES galaxies via the completion of the IRAS 60^tm luminosity function with pure luminosity 
evolution up to z = 2 (§3.3). The (pink) dashed-line shows the SFRD inferred from 24/jm-selected luminous and ultraluminous infrared 
galaxies (Le Floc'h et al. 2005), which follow a similar luminosity evolution to that presented here for the SHADES galaxies. The (purple) 
thick dash-dotted line shows the SFRD for the sample of SCUBA galaxies from Chapman et al. (2005). The thin dash-dotted line, which 
is a factor of ~ 3 higher, is an estimation of the contribution of 850/^m-selected galaxies down to the ~ 1 mjy level. The SFRDs for 
optical/UV-selected starbursts are shown as small triangles, and are taken from Lilly et al. (1996 - upward-pointing triangles), Connolly 
et al. (1997 - downward-pointing), Steidel et al. (1999 - left-pointing) and Giavalisco et al. (2004 - right-pointing). The optical/UV 
data are shown with and without corrections for dust-obscuration as empty and solid triangles respectively, and all data have been 
homogenized to the same set of parameters and corrected to complete a Schechter luminosity- function (Giavalisco et al. 2004). All 
estimates have been converted to the same SFR/Lir factor and cosmological model described in § 1. 



follows (1 + z) 3 at 2 ^ 2, and then maintains a constant 
level for z > 2. This evolutionary form provides an ad- 
equate description of the 850pm and 1.1mm counts (e.g. 
Scott et al. 2002, Greve et al. 2004) and is also supported 
up to z — 1 by the evolution of the luminosity function 
of 24/im-selected galaxies (Le Floc'h et al. 2005). By imple- 
menting this luminosity-function correction the SFR density 
increases up to a maximum factor of 2 at z > 2. 

The contribution of SHADES galaxies to the global 
star formation rate density of the Universe is comparable to 
the contribution of starbursts selected at optical/UV wave- 
lengths at 1 <, z <, 4 before the latter are corrected for dust 



extinction. It is important to recall that the Lfir/SFR fac- 
tor we have adopted could be in error by ~ ±30 per cent, 
and that SED differences could also account for an increase 
of ~ 40 per cent. These uncertainties have not been carried 
into the estimation of error-bars in figure 9. If we complete 
the luminosity function of SHADES galaxies towards lower 
luminosities, the FIR star formation rate traced by ultra- 
luminous starbursts is still a factor of 1.2 to 2 lower than 
that of optical/UV starbursts that have been corrected for 
intrinsic dust extinction. 

The recent demonstration that the contribution of lu- 
minous and ultraluminous IR galaxies dominates the SFR 
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density at z ^ 1 (Le Floc'h et al. 2005) suggests that even 
if we correct for the incomplete sampling of the sub-mm 
galaxy luminosity function, a bright SHADES sub-mm sur- 
vey could be missing the integrated contribution of dusty 
starbursts to the global SFR by a large factor (as high as 
~ 7), and thus dusty starbursts could indeed prove to be a 
significant mode of the star formation of the Universe (Blain 
et al. 1999). 

The contribution of the fainter (< 3mJy) 850/im sub- 
mm galaxy population to the star formation history of the 
Universe at z > 1 remains unconstrained at present, since 
detecting faint sub-mm galaxies has been restricted to a few 
strongly-lensed fields (e.g. Smail, Ivison & Blain 1997) and 
extremely-deep confusion-limited pencil observations (e.g. 
Hughes et al. 1998). Our estimations presented in figure 9 
should therefore be considered an educated estimate of how 
ultraluminous IR-submm galaxies trace the star formation 
history of the Universe. 

Accurate measurements of the surface density and red- 
shift distribution of the entire sub-mm galaxy population 
(and the faintest galaxies in particular) that contribute the 
complete sub-mm to FIR extragalactic background require 
deeper and larger mm and sub-mm surveys than are cur- 
rently possible. The anticipated continuum and spectro- 
scopic surveys with SCUBA-2 (Holland et al. 2006), the 
Large Millimetre Telescope (LMT, Serrano et al. 2006), 
and the Atacama Large Millimeter Array (ALMA, Beasley, 
Murowinski & Tarenghi 2006), for example, will provide 
suitable data. In the meantime, however, it is still possi- 
ble that the redshift distribution of the more populous and 
fainter (possibly extremely high-redshift or alternatively less 
luminous) sub-mm galaxies is significantly different to those 
galaxies identified in the SHADES survey. 



4 SUMMARY OF CONCLUSIONS 

• We have derived the photometric redshift distribution 
of SHADES sources with de-boosted 850/im flux densities 
> 3 mjy towards the Lockman Hole (LH) and the Sub- 
aru/XMM Newton Deep Field (SXDF) using rest-frame ra- 
dio to FIR photometry. The redshift distribution of the 
radio-detected sub-mm sources peaks at z ~ 2.4 with 50 per 
cent of the population between redshifts 1.8 and 3.1. 

• The combined redshift distribution of SHADES sources 
with robust radio counterparts, ~ 60 per cent of the pop- 
ulation, has a distribution which is qualitatively consistent 
with the distribution of rest-frame UV-optical spectroscopic 
redshifts published by Chapman et al. (2005). 

• We find a small [5z ~ 0.5), but significant, difference 
between the peaks of the photometric-redshift distributions 
of the LH and SXDF, which can be attributed to differences 
in the sensitivities of their respective radio maps. Intrinsic 
field-to-field redshift variance is also expected, and is charac- 
terized according to a variety of models (8z w 0.25 — 0.55). 
This drives us to the conclusion that the incomplete area 
(~ 720 sq. arcmins) observed by SCUBA, despite being the 
largest sub-mm survey to date, may still be too small to be 
a representative sample of the bright sub-mm galaxy popu- 
lation. 

• The complete redshift distribution of all SHADES 
sources, including those sub-mm sources without detections 



at radio wavelengths (for which we adopt a variety of pos- 
sibilities that describe their unconstrained distribution of 
redshifts) still maintains the peak (mode) of the bright sub- 
mm galaxy redshift distribution at ~ 2.4. We have consid- 
ered a variety of priors that describe the unconstrained red- 
shift distributions of the sub-mm sources without radio de- 
tections. In the most extreme solutions, distributing these 
sources with equal probability between their lower redshift- 
limit and z = 2 or z — 5, the bulk of the sub-mm population 
(50 per cent interquartile) lies in the range 1.6 < z < 2.6 or 
2.1 < z < 3.4, respectively. 

• The combined SHADES LH and SXDF redshift 
probability-density distribution is compatible, within the 
uncertainties of our analysis, with the semi-analytical model 
for the joint formation of spheroids and QSOs of Granato 
et al. (2004) and Silva et al. (2005). If sources detected only 
at 850/im are also introduced into the redshift probability- 
density, with other priors than those illustrated here, then 
the hydrodynamical model of Muanwong et al. (2002) and 
phenomenological model of van Kampen (2004) and van 
Kampen et al. (2005) could also be in agreement with the 
observations. These compatible models, which are physically 
quite distinct, predict different clustering properties for the 
SHADES galaxies that could allow further discrimination 
between them (van Kampen et al. 2005). A detailed study 
of the clustering properties of SHADES galaxies will be the 
topic of a further paper. 

• The bright SHADES galaxies contribute to the SFR 
density of the Universe with ~ 0.01 to 0.03 M yr _1 Mpc~ 3 
in the redshift interval 1 <, z <, 5, and reach the levels of 
the dust-uncorrected Lyman Break Galaxy population (Gi- 
avalisco et al. 2004). The SFR density of dust-enshrouded 
starburst galaxies traced by ultraluminous SHADES galax- 
ies, and completing the luminosity function to lower lu- 
minosity galaxies, is estimated to be a factor of 2 larger. 
This is still a factor of 1.2 to 2 lower than the optical/UV- 
selected starburst galaxy samples that include the lat- 
est dust-correction estimates. The current SHADES sur- 
vey and complementary multiwavelength data, however, 
cannot characterize the bulk of the rest-frame FIR emis- 
sion arising from these lower luminosity galaxies. A more 
statistically-complete measurement of the universal history 
of star-formation from powerful dusty, optically-obscured 
galaxies awaits the commissioning of future large-aperture 
single-dish and interferometric submillimetre and millimetre 
telescopes targetting suitable extragalactic fields that have 
the necessary multi-wavelength ancillary data. 
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Table 1. Photometric redshifts for SHADES sources in the LH field based on the 1.4GHz/850/xm spectral index. The columns give: 
(1) name of the source; (2) z photometric redshift using the prescription of Carilli & Yun (1999, 2000); (3) 2 ^of photometric 
redshift using the template collection of Aretxaga et al. (2003, 2005), the 90% confidence interval is given in parenthesis; (4) notes 
on which radio counterpart (from Ivison et al. 2007) is used in the computation of photo-z, in case of ambiguity (N for Northern 
component, S for Southern component, etc, or 'coadded' if the flux densities from all components are summed); and (5) z S pec, 
spectroscopic redshift taken from the literature, where sources for which the redshifts are in parenthesis have reported ambiguities 
in their radio/optical counterpart associations, or where the redshifts are otherwise under scrutiny. The references for the spec-z 
(as a superscript of the values) and any debate about them (after the parenthesis, where it applies) are as follows: 1.- Chapman et 
al. 2005; 2.- Ivison et al. 2005; 3.- Greve et al. 2005; 4.- Chapman et al. 2003; 5.- Swinbank et al. 2005; 6.- Chapman et al. 2002; 
7.- Almaini et al. 2005; 8.- Kovacs et al. 2006. 
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a This source has a robust 1.4GHz association in the dataset of Ivison et al (2002), but it is below the robustness level adopted for 
the analysis in this paper, and thus we will make use of the revised 1.4GHz photometry of Ivison et al. (2007) as an upper limit. 
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Table 1. (cont.) 
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Table 2. Same as table 1 for SHADES sources in the SXDF field. 
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5 


9) 


SHADES 


J021756-045806 


(SXDF850 


9) 


^ 2.1 




1.6 






SHADES 


J021825-045557 


(SXDF850 


10) 


2 3+ 1 ' 3 


1 6+ 1 ' 4 


(1.0- 


-4 


8) 


SHADES 


J021725-045937 


(SXDF850 


11) 


2.8±l;S 


2 1+ 1 - 6 


(1.0- 


-5 


6) 


SHADES 


J021759-050503 


(SXDF850 


12) 


3.6±?-§ 




(1.5- 


6 


7) 


SHADES 


J021819-050244 


(SXDF850 


14) 


3.4±i: D 


2-6±i : g 


(1.2- 


6 


5) 


SHADES 


J021815-045405 


(SXDF850 


15) 


^ 2.5 




2.0 






SHADES 


J021813-045741 


(SXDF850 


16) 


3 5-I- 2 - 1 


2 6+ 1 - 9 


(1.5-6 


5) 


SHADES 


J021754-045302 


(SXDF850 


17) 


2.7 




2.2 






SHADES 


J021757-050029 


(SXDF850 


18) 


3.9±2;5 


2-6±g;g 


(1.5- 


6 


5) 


SHADES 


J021828-045839 


(SXDF850 


19) 


2 2+ 1 ' 2 




(0.8- 


-4 


6) 


SHADES 


J021744-050216 


(SXDF850 


20) 


^ 1.7 




1.4 






SHADES 


J021742-050427 


(SXDF850 


21) 


0.9±0;f 


o-6±8:§ 


(0.0- 


-2 


2) 


SHADES 


J021800-050741 


(SXDF850 


22) 


1.7 




1.8 






SHADES 


J021742-050545 


(SXDF850 


23) 




2 1 + 1 ' 4 


(1.0- 


-5 


°) 


SHADES 


J021734-050437 


(SXDF850 


24) 


3.4± 2 ° 


2 l-l- 1 ' 8 
2 6+ 2 ' 8 


(1.0- 


6 


1) 










3-7±?1 


(1.5- 


-7 


0) 


SHADES 


J021812-050555 


(SXDF850 


25) 


^ 1.6 




1.0 






SHADES 


J021807-050148 


(SXDF850 


27) 


i k_i_0.8 
i - O=t 0.B 


1 l+°' 9 

1 6+ 1 ' 4 
1 - D=C 0.6 


(0.5- 


-3 


3) 


SHADES 


J021807-045915 


(SXDF850 


28) 


2 3+ 1 ' 3 


(0.5- 


-4 


5) 










1 7+ 10 


i - ±:c 0.3 


(0.2- 


-3 


5) 


SHADES 


J021816-045511 


(SXDF850 


29) 


1 6± ' 9 

0=, -Q.6 


1 1±°' 9 


(0.6- 


-3 


5) 


SHADES 


J021740-050116 


(SXDF850 


30) 


4.4±|;J 


3.1±J;I 


(1.8- 


-7 


7) 


SHADES 


J021736-045557 


(SXDF850 


31) 


3-2±i:§ 


2-i±i:i 


(1.5- 


6 


2) 


SHADES 


J021722-050038 


(SXDF850 


32) 


2.1 




1.5 






SHADES 


J021800-045311 


(SXDF850 


35) 


3.3±?;g 


2 1+ 21 


(1.0- 


6 


1) 


SHADES 


J021832-045947 


(SXDF850 


36) 


^ 1.9 




1.8 






SHADES 


J021724-045839 


(SXDF850 


37) 


3.2±?;g 


2 1 + 2 ' 1 


(0.8-6 


1) 


SHADES 


J021825-045714 


(SXDF850 


38) 


2 7+ 1 ' 6 




(0.5- 


5 


4) 


SHADES 


J021750-045540 


(SXDF850 


39) 


1.7 




1.5 






SHADES 


J021729-050059 


(SXDF850 


40) 


2.9±;-« 


2 l+ 2 '° 


(1.0- 


-5 


9) 


SHADES 


J02 1829-050540 


(SXDF850 


45) 


^ 4.4 




3.3 






SHADES 


J021733-045857 


(SXDF850 


47) 


1 5+ - 9 
1 - OzSz 0.6 


1 1+°- 8 
i - 1± 0.6 


(0.2- 


-3 


1) 










2.6±if 


2 1 + 1 ' 4 


(0.6- 


-5 


2) 










1 3+° :7 


1 1+°' 7 


(0.2- 


2 


9) 










2 2+ 1 ' 2 


1 6+ 1 ' 3 


(0.5- 


4 


4) 


SHADES 


J021724-045717 


(SXDF850 


48) 


> 2.6 


> 


1.8 






SHADES 


J021820-045648 


(SXDF850 


49) 


^ 1.3 




1.0 






SHADES 


J02 1802-045645 


(SXDF850 


50) 


3.6±?;3 


2-l±^ 


(1.2- 


-6 


9) 


SHADES 


J02 1804-050453 


(SXDF850 


52) 


2-o±Ji 




(0.5- 


4 


1) 










i-5±R:S 


1 1±°- 9 

2 1+ 2 ' 1 


(0.2- 


-3 


1) 


SHADES 


J021752-050446 


(SXDF850 


55) 


3-0±L 7 


(0.8- 


6 


1) 


SHADES 


J021750-050631 


(SXDF850 


56) 


1.2 




0.6 






SHADES 


J021745-045750 


(SXDF850 


63) 


^ 1.6 




1.5 






SHADES 


J021807-050403 


(SXDF850 


65) 


^ 2.2 




1.2 






SHADES 


J021751-050250 


(SXDF850 


69) 


1.5 




1.0 






SHADES 


J021811-050247 


(SXDF850 


70) 


^ 1.9 




1.2 






SHADES 


J021821-045903 


(SXDF850 


71) 


^ 1.7 




1.2 






SHADES 


J021758-045428 


(SXDF850 


74) 


2-9±i;I 


2.1±J;f 


(0.8- 


-5 


8) 


SHADES 


J021755-050621 


(SXDF850 


76) 


1.9 




1.5 






SHADES 


J021736-050432 


(SXDF850 


77) 


2-9±i; s 7 


2 1+ 2 - 1 


(1.0- 


-6 


1) 


SHADES 


J021817-050404 


(SXDF850 


86) 


^ 1.5 




1.0 






SHADES 


J021800-050448 


(SXDF850 


88) 


^ 1.0 




1.0 






SHADES 


J021734-045723 


(SXDF850 


91) 


1.5 




1.2 






SHADES 


J021733-045813 


(SXDF850 


93) 


^ 1.9 




0.8 







NW 
N 



N 
S 



N 

N+S 



NE 
SE 

NE+SE 
W 



E 
ES 
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Table 2. (cont.) 



object 


2 CY 
phot 


Z phot n ° teS Z *P°c 


SHADES J021740-045817 (SXDF850.94) 
SHADES J021741-045833 (SXDF850.95) 
SHADES J021800-050212 (SXDF850.96) 
SHADES J021756-045255 (SXDF850.119) 


^ 1.5 
^ 1.6 
3.5±?;i 
3-4±?| 


> 1.2 

> 1.0 
3.1±H (1.0-6.3) 
2.6±} : J (0.8-6.2) 



Table 3. Photometric redshifts for the sources with additional radio-mm-FIR data reported in the literature from other surveys. 
The photometry has been complemented at other wavelengths, while preserving the SHADES 450/850/jm (Coppin et al. 2006) and 
1.4GHz (Ivison et al. 2007) determined flux densities. The first column gives the source-name; the second column gives the most 
probable mode and error bars based on the 68% confidence interval of the mode calculation (in parenthesis the 90% confidence 
interval); the third and fourth columns respectively give the bands in which the source is detected at a ^ 3cr level, and at which 
upper-limits are used for the computation of the photomctric-rcdshifts; the fifth column provides the published references (and 
alternative sources-names in overlapping surveys) to the additional photometry, using the following syntax: LE850.X from Scott et 
al. 2002, 1100.x from Laurent et al. 2005, 2006 and 1200.x from Greve et al. 2004 and Ivison et al. 2005). 



object 


■^phot 




J? 3a detections 


< 3cr 1 upper limits 


notes 


Lock850.1 




(2.2- 


-3.8) 


350,850Aim,l.l,1.2mm,1.4GHz 


175,450/mi,3.3mm,5GHz 


LE850.1, 1200.5, 1100.14 


Lock850.2 


2-9±8:f 


(2.5- 


-3.8) 


350,850/im, 1.1,1.2mm, 1.4GHz 


450/xm 


LE1100.1, 1200.4, SW-1.4GHz 




2-9±8:I 


(2.8- 


-3.9) 






NW-1.4GHz 


Lock850.3 




(2.5- 


-4.2) 


350,850/im, 1.1,1.2mm, 1.4GHz 


450/zm, 5GHz 


LE850.2, 1100.8, 1200.1, S-1.4GHz 




2.6±0;3 


(2.5- 


-3.8) 






coadded-1.4GHz 


Lock850.4 


1.6±g;? 


(1.5- 


-4.8) 


850^tm, 1.2mm, 1.4GHz 


450/zm, 5GHz 


LE850.14, 1200.3 


Lock850.12 




(2.2- 


-3.0) 


350,850/im, 1.1,1.2mm, 1.4GHz 


450/im, 5GHz 


LE850.16, 1100.16, 1200.6 


Lock850.14 


2.6±8;? 


(2.2- 


-3.7) 


350,850£tm, 1.1,1.2mm 


450Atm, 1.4,5GHz 


LE850.6, 1100.5, 1200.10^. 


Lock850.16 


i-9±8 :i 


(1.5- 


-3.2) 


850Atm, 1.2mm,1.4GHz 


450/mi, 5GHz 


LE850.7 


Lock850.17 


2 5+ - 6 


(2.0- 


-5.9) 


850jUm, 1.2mm, 5,1.4GHz 


450/im 


LE850.3, 1200.11 


Lock850.18 


3 1+2-9 


(2.3- 


-6.0) 


850/im, 1.2mm, 1.4GHz 


450/im 


LE1200.9 


Lock850.27 


4 6+ 1 - 4 


(4.0- 


-6.0) 


850/mi, 1.1,1.2mm, 1.4GHz 


450/mi, 5GHz 


LE1100.4, 1200.7 


Lock850.33 


3 6+ - 7 


(2.4- 


-4.8) 


850^m, 1.2mm, 1.4GHz 


450/xm 


LE850.18, 1200.12 


Lock850.41 




(3.2- 


-4.4) 


350,850/im, 1.1,1.2mm, 1.4GHz 


450/xm 


LE850.8, 1100.17, 1200.14 


Lock850.76 


4.6±l;f 


(3.0- 


-6.0) 


850^m,l.lmm,1.4GHz 


450/im 


LE1100.15 



a This source has a robust 1.4GHz association in the dataset of Ivison et al (2002), but it is below the robustness level adopted 
for the analysis in this paper, and thus we will make use of the 1.4GHz photometry of Ivison et al. (2007) as an upper limit. 
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Table 4. Photometric redshifts for SHADES sources in the LH field based on the 850/xm and 1.4GHz data and 450/xm upper limits 
determined by SHADES. The columns give: (1) name of the source; (2) z S ^ ot survival analysis solution; (3) z^^. Monte Carlo 
solution; and (4) notes on associations. 



object 



Z SA 
phot 



Z MC 
phot 



notes 



Lock850.1 
Lock850.2 

Lock850.3 

Lock850.4 

Lock850.5 

Lock850.6 

Lock850.7 

Lock850.8 

Lock850.9 

Lock850.10 

Lock850.11 

Lock850.12 

Lock850.13 

Lock850.14 

Lock850.15 

Lock850.16 
Lock850.17 
Lock850.18 
Lock850.19 
Lock850.21 
Lock850.22 
Lock850.23 
Lock850.24 
Lock850.26 
Lock850.27 
Lock850.28 
Lock850.29 
Lock850.30 
Lock850.31 
Lock850.33 
Lock850.34 
Lock850.35 
Lock850.36 
Lock850.37 

Lock850.38 
Lock850.39 
Lock850.40 
Lock850.41 

Lock850.43 
Lock850.47 
Lock850.48 
Lock850.52 
Lock850.53 
Lock850.60 
Lock850.63 
Lock850.64 
Lock850.66 
Lock850.67 
Lock850.70 
Lock850.71 



2.1±g;5 (1.5-5.9) 
3.1±g : f (2.4-6.8) 
3.6±g : g (2.5-6.8) 
4-l±gJ (2.8-7.7) 
3.1±J;g (2.0-6.1) 
2.1±J-| (1.5-5.8) 

> 2.9 
4.1±g-| (3.0-8.0) 
3.1±j| (2.0-6.3) 

^ 2.5 
2.1±g-; (1.5-5.8) 
3.4±§ : § (2.2-7.0) 

> 2.2 
2.6±J;| (1.5-6.2) 

^ 1.5 

^ 2.2 
2.6±i-| (1.2-5.9) 
3.1±2 ° (1.8-6.8) 
1.6±i| (1.0-4.7) 
1.6±i-2 (1.0-4.2) 
3.1±2-i (1.5-6.1) 

> 1.8 
^ 1.5 
^ 2.0 

> 1.6 
2.6±J;f (1.2-5.8) 
3.1±J : f (1.5-7.2) 
3.9±^ (1.8-6.5) 

^ 2.0 

> 2.2 
l.l±°'f (0.5-3.2) 
2.6±J-jj (1.5-6.4) 
2.1±J-| (1.2-5.4) 
3.4±J : g (2.0-6.5) 

^ 2.0 

^ 2.5 
2.9±J 
3.9±? 
2.6±J;f (0.1 

> 2.2 
2.6±g;| (1.2-6.3) 



1.2.6 

.1 



-6.2) 



2 1+ 1 - 4 
1 4± 1 ' 4 
3.6±P 



(1.2-5.3) 
(1.0-4.7) 
(1.5-6.6) 
$ 1.2 
(0.5-5.7) 
(0.5-5.5) 
^ 1.5 
^ 1.2 
(1.5-6.4) 

> 1.6 
^ 1.5 
^ 1.0 

> 1.0 
1.6±J;| (0.5-4.3) 



2-i±i:f 



)_2.1 
c 0.6 



(2.6- 



3.9±g;$ (3.2- 



2-9±8:f 

3 R+0.4 

uO.fi 
c 0.4 
2 6+ ' 4 

2-i±8:? 

3.6±l;0 

2-a±8:l 



(1.8- 
3.0 



2.4± 



l0.3 
-0.4 

3 1+ 09 
d - i=t 0.3 

2 fi+ ' 4 

2.4±° 
2.9±° 

3-l±g 
2.9±° 
2.9±« 



i 

7 
3 
4 
4 
5 

2 4+ - 4 
3.6±8;J 
3.9±if 

2 1± 01 
2 6+ d ' 1 

2 1+ - 7 

3 1+°- 6 

4.5± 1 - 
4.6±|i 

i 0.6 
c 0.2 
3 fi+0- 4 

2 q+°' 5 

2 4+ 8 ' 8 

2 4+0.5 
2 1+ - 1 

2 6+ - 4 

Z.DICq 4 

2 9+ 01 



2.5) 




3.8) 


SW 


4.2) 


NW 


4.8) 


S 


3.2) 


coadded 


3.1) 


coadded 


4.8) 




3.4) 




3.2) 




4.7) 




3.8) 





2.6±^ 



2.5 



2.5 
(2.2- 
1.5 
2.2 

(2.0-3.2) 

(2.2-3.8) 

(2.0-3.4) 

(2.3-3.2) 

(2.2-3.7) 

1.5 

1.0 

2.0 

1.5 

(2.0-3.1) 
(2.2-3.9) 
(1.8-6.5) 
2.0 
2.1 
(1.8- 
(2.0- 
(1.8-3.2) 
(2.6-3.8) 
2.0 
2.5 

(4.3-5.8) 
(4.5-5.8) 
(2.0-2.3) 
2.0 

(2.0-3.2) 
(2.5-4.0) 
(2.3-3.5) 
(2.2-3.8) 
1.5 

(2.1-3.0) 
(1.9-2.2) 
1.5 
1.5 

(2.1-3.2) 

1.5 

1.5 

1.0 

1.0 

(2.0-3.2) 



2.5) 
3.1) 



no 1.4GHz^ 

coadded 

S 



N (P = 0.013) 

S (adopted P = 0.078) 



S 

N+S 

adopted P = 0.060 
adopted P = 0.068 



a This source has a robust 1.4GHz association in the dataset of Ivison et al (2002), but it is below the robustness level adopted 
for the analysis in this paper, and thus we will make use of the 1.4GHz photometry of Ivison et al. (2007) as an upper limit. 
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Table 4. (cont.) 



object 



Z SA 
phot 



Z MC 
phot 



notes 



Lock850.73 


2 1+ 2 ' 6 


(1.0-6.2) 


2-4±;i: 2 


(2.0-2.5) 


N 




2 1+ 1 - 4 


(0.7-5.2) 


2-i±8:i 


(1.9-2.5) 


coadded 


IjOCKoOU- I O 




1 1 


-> 






Lock850.76 


2-i±5:i 


(0.8-5.5) 


2-i±8:5 


(2.0-2.5) 




Lock850.77 


2.9±i;[ 


(2.8-6.9) 


2.6±g;S 


(2.2-3.8) 


S 




1 5+ 1 - 6 
l.Oig 4 


(0.8-4.8) 


2.8±°;t 


(2.0-3.2) 


N+S 


Lock850.78 




1.1 




1.5 




Lock850.79 


2-6±§;| 


(1.2-6.3) 


2.4±8;S 


(2.0-3.2) 


adopted P = 0.064 


Lock850.81 




1.9 




2.0 




Lock850.83 




0.8 




1.5 




Lock850.87 


1 6+ 10 


(0.8-4.0) 


2.4±o : ; 


(1.9-3.0) 




Lock850.100 




2.1 




2.0 
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Table 5. Photometric redshifts for SHADES sources in the SXDF field based on the 850/im and 1.4GHz data and 450/im upper 
limits, and when significant, complemented with Spitzcr photometry at 70 and 160/im. Columns arc as in table 4. 



object 

SXDF850.1 
SXDF850.2 
SXDF850.3 
SXDF850.4 
SXDF850.5 
SXDF850.6 

SXDF850.7 

SXDF850.8 

SXDF850.9 

SXDF850.10 

SXDF850.il 

SXDF850.12 

SXDF850.14 

SXDF850.15 

SXDF850.16 

SXDF850.17 

SXDF850.18 

SXDF850.19 

SXDF850.20 

SXDF850.21 

SXDF850.22 

SXDF850.23 

SXDF850.24 

SXDF850.25 
SXDF850.27 
SXDF850.28 

SXDF850.29 
SXDF850.30 
SXDF850.31 
SXDF850.32 
SXDF850.35 
SXDF850.36 
SXDF850.37 
SXDF850.38 
SXDF850.39 
SXDF850.40 
SXDF850.45 
SXDF850.47 



SXDF850.48 
SXDF850.49 
SXDF850.50 
SXDF850.52 

SXDF850.55 
SXDF850.56 
SXDF850.63 
SXDF850.65 
SXDF850.69 
SXDF850.70 
SXDF850.71 
SXDF850.74 
SXDF850.76 
SXDF850.77 
SXDF850.86 



'phot 



'phot 



notes 



2 9+ 2 ' 
z - a *0.6 

2.4±°;f 

9 4+1-5 
z - 4 *0.6 
1 1+1-° 

1 4± ' 6 

2 4+1- 9 
z - 4 *0.6 

2 9+1- 6 
z - a *0.9 

2 4+ 2 - 1 

2 4+1- 6 
z - 4 *0.6 

2-4±J:S 



2 4+1- 9 
z - 4 *0.6 



2 4+ 2 - 
z - 4 *0.6 

1-6±H 



0.6±°; 9 



1 9+1- 4 
!- a *0.6 



2.4+i 



_1.6 
-0.9 
2.9±l;f 

1 2+11 
l- z *0.4 

1 l+°' 9 
J - 1=c g.4 

2.9±?;f 

2 1 +1- 6 



9 4+1-6 
z - 4 *0.9 



2 

''■ 1=t 0.5 

1 9+ 1 ' 8 
!- a *0.6 



2 1+ 1 ' 9 
z -l*0.6 



1 1+°' 8 

2 1+ 1 - 4 

1 1+ - 7 

1-1*0.6 

1 6+ 1 - 3 
!- D *0.6 

2.4± 2 -0 
1.6±, 



.1 

-0.6 

1 l+°- 9 
1-1*0.6 

2.4±;-f 
2.i±;-f 

2 1+ 1 - 
z -l*0.6 



1.8-6.3 

1.0-4. 

1.5-5.9 

0.5-3.4 

0.5-2.8 

1.3-6.0 

1.5-6.0 

1.8-6.6 

1.2-5.8 

.8 

0.8-4.6 
1.0-5.6 
1.2-6.3 
1.2-6.3 
.8 

1.2-6.2 
.0 

1.5-6.1 
0.8-4.6 
4 

0.0-2.2 
.8 

1.0- 5.0 

1.1- 6.0 
1.5-6.8 
.0 

0.5-3.3 
0.5-4.4 
0.2-3.2 
0.6-3.5 

1.2- 6.8 
1.5-6.1 
.5 

1.2-6.2 
.8 

0.8-5.9 
0.3-5.0 
.2 

1.0-5.8 
1.8 

0.2-3.1 
0.6-5.2 
0.2-2.9 
0.5-4.4 
.8 
'.5 
1.2-6.6 
0.5-4.1 
0.2-3.1 
0.5-5.6 

0.6 

1.0 

1.2 

0.8 

1.0 

0.8 

(0.9-5.8) 
1.2 

(1.0-6.0) 
1.0 



2-6±8:t 

2 1+0- 3 
2 1+ - 4 

1 4+ - 4 

2 4+ 1 - 2 
2 6+ - 4 

Z.D=C 4 

9 4+0.5 

^■^o.i 

2 6+ 1 - 3 

2 6+ - 3 

2 4+ - 4 
fi- 4 
9 4+0.3 

2 4+ - 4 

2-4±8:S 
2-9±g:I 

2 4+ - 4 
0.5±°; 4 

9 4_|_0.9 
9 4+0.6 

2 9+0- 1 
Z - y=t 0.6 



3.6±° 
1.4±« 
1.1±8 
2.1±° 
2.9±0 
2.6±° 



5 
4 

9 4+0-6 



2 1+ - 1 
1 9+0- 2 



2-4±8:l 
i.9±8:? 

i-e±8:f 

2 1+0- 1 

2 9+0 1 

2 1+0-1 

2 1+ - 1 

2 1+0- 1 
^.rzn n 



1± ().l 



1+0.1 
• ±=c ().l 



(2.2-3.4 
(1.8-2.9 
(2.0-2.6 
(1.8-2.7 
(0.6-2.0 
(2.1-3.8 
(2.2-3.7 
(2.2-3.4 
(2.3-4.0 
1.5 

(1.9-3.2 
(2.0-3.4 
(2.0-3.0 
(2.0-3.1 
2.0 

(2.0-3.2 
2.0 

(2.2-3.7 
(2.0-3.2 
1.5 

(0.0-1.2 
2.0 

(2.0-3.5 
(2.2-3.7 
(2.2-3.7 
1.0 

(2.4-3.8 
(0.6-2.2 
(0.5-1.5 
(1.8-2.4 
(2.2-3.8 
(2.2-3.7 
1.5 
(2.1 
2.0 
(2.0-2.4 
(1.7-2.2 
1.5 

(2.0-3.2 
3.0 

(1.6-2.2 

(1.9-2.2 

(1.5-2.1 

(1.8-2.4 

2.0 

1.0 

(2.2-3.6 

(1.6-2.2 

(1.5-2.2 

(1.8-2.2 

1.0 

1.0 

1.5 

1.0 

1.0 

1.0 

(1.8-2.2) 
1.5 

(1.8-2.2) 
1.0 



-3.2 



70,160/im included 

NW 

N 



70,160/im included 



N, 70,160/im included 
N+S 



NE 
SE 

NE+SE 
W 



E 
ES 



70, 160/im included 
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Table 5. (cont.) 



object 








notes 


SXDF850.88 


^ 1.0 


^ 1.0 






SXDF850.91 


^ 1.0 


^ 1.5 






SXDF850.93 


^ 0.8 


^ 0.5 






SXDF850.94 


^ 1.2 


^ 1.5 






SXDF850.95 


^ 1.0 


^ 1.0 






SXDF850.96 


2.4±J f (1.0-6.1) 


2.4±g f (2.2- 


-3.5) 




SXDF850.119 


2.2±i | (0.0-4.5) 


1.9±{& (1.2- 


-2.2) 


70,160/^m included 
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